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^J^ ' We give an example of modeling phcnomenological heavy-quark potentials in a five- 

f ^ . dimensional framework nowadays known as AdS/QCD. In particular we emphasize the ab- 

"^ ' sence of infrared renormalons. 
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\r^ • 1 Introduction 

C^ ' Heavy-quark potential is one of the basic observables relevant to confinement. It has been 

measured in great detail in lattice simulations^ and the results reveal a remarkable agreement 
with the so called Cornell potential |U 

V{r) = -- + -^+C, (1) 

where the coefficients are adjusted to fit the charmonium spectrum 

K«0.48, asa2.34GeV-^ C = -0.25GeV. (2) 

As follows from above, n and -^ can be interpreted as gCts and the string tension, respectively. 
One of the implications of the AdS/CFT correspondence ^3, is that it resumed interest in 
finding a string description of strong interactions. For the case of interest, let us briefiy mention 
a couple of results. 
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First, in the approach called usually gauge/string duality one tries to keep the underlying 
string structure. As a consequence, the theory is ten dimensional and its reductions to five dimen- 
sions in general contain additional higher derivative terms (stringy a' corrections).^ According 
to 131 , the expectation value of the Wilson loop is given by 

(VF(C))~e-^ (3) 

where S is an area of a string world-sheet bounded by a curve C at the boundary of AdS space. ^ 
Second, in more phenomenological approach called AdS/QCD one starts from a five dimen- 
sional effective field theory somehow motivated by string theory and tries to fit it to QCD as 
much as possible. It was recently pointed out [Sl[7j that asymptotic linearity of Regge trajectories 
arises for some backgrounds. Such backgrounds reduce to the standard AdS background in the 
UV but differ from it in the IR. The latter turns out to be crucial for linearity. In this case, it 
would be natural to expect that the interquark interaction would include the dominant Coulomb 
term at short distances as well as the dominant linear term at large distances. 

In this note we explore this expectation in the context of (jSJ. For simplicity we will concen- 
trate here on the case of [Jj. However, it should not be hard to adapt the arguments to 6 . It is 
worth noting a recent attempt to derive the Cornell type potential within the model based on a 
truncated AdS space |^.^ 

Before proceeding to the detailed analysis, let us set the basic framework. We consider the 
following Euclidean background metric 

ds"^ = GnrndX-^dX"" = R^ ^ idx'dx' + dz^) , h = Q^^'' , (4) 

where i = 0, . . . , 3. In the region of small z the metric behaves asymptotically as Euclidean AdSs, 
as expected. We also take a constant dilaton. Note that the use of the Euclidean signature for 
the background metric slightly modifies h. So, the exponent has an opposite sign to that of |7|. 

2 Calculating the Potential 

Given the background metric, we can attempt to calculate the corresponding potential. In doing 
so, we adapt the conjecture Q to AdS/QCD. 

We consider a rectangular Wilson loop C living on the boundary [z = 0) of five dimensional 
space as shown in Fig.l. The quark and antiquark are set at x = | and a; = — |, respectively. 
As known, taking the limit T ^f oo allows one to read off the energy of such a pair from the 
expectation value of the Wilson loop namely, {W{C)) ~ e"^-^^''). 

Now we are ready to evaluate the expectation value of the loop. To this end, we make use of 
the Nambu-Goto action equip with the background metric @ 



1 ' .2. 
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d'i Wdet G^^do^X^dpX"^ . (5) 



^These corrections are of order —^- Thus, they might be relevant at Nc = 3. 

^The literature on the Wilson loops within the AdS/CFT correspondence is very vast. For a discussion of this 
issue, see, e.g., and references therein. 

^In this case there is a subtle point. The use of two different solutions leads to a discontinuity in the interquark 
force. 



Figure 1: A loop C. 



Next, we choose £,^ =t and ^^ = x. This yields 



27r _L z^ 



i?2 



where g = — . A prime denotes a derivative with respect to x. 
Now it is easy to find the equation of motion for z 

zz" + (2 - cz") (1 + {z'f) = 

as well as the first integral 



h 



const . 



(6) 



(7) 



(8) 



z^^l + {z>f 

The integration constant can be expressed via the maximum value of z. On symmetry grounds, 
z reaches it at x = 0. By virtue of (ISJ, the integral over [—§,§] of dx is equal to 



r = 2^ f' dv ^2e^^(i-^') f 1 - v'e^^'-"'A ^ 



(9) 



where v = j-, X = czq, and zq = z\x=o- 

At this point a comment is in order. A simple analysis shows that the integral © is real for 
A < 2. It develops a logarithmic singularity at A = 2 and becomes complex for larger A. Hence, 
there exists the upper bound on the maximum value of z 



zo < 



(10) 



Note that in the limit as c goes to zero zq is not bounded, as should be for the AdS space. As 
function of z, the effective string tension reaches its minimum at z = zq. Thus, there exists a 
kind of horizon which is a generic feature of confining theories. 

Now, as in i;4, we will compute the energy of the configuration. First, we reduce the integral 
over x in Eq.® to that over z. This is done by using the first integral ((HI). Since the integral 
is divergent at z = due to the factor z~'^ in the metric Q, in the process we regularize it by 



imposing a cutoff e. Then we replace z with w as in (jU)). Finally, the regularized expression takes 
the form 

E^ = ^,l^- [\vv-'e'2^^'(l-v'e'('-^'^y\ (11) 

Its e-expansion is simply 
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E^ = -^ + E + Oie), (12) 
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where 




E = ^J- \ -1+ I dvy-^ 



dvv e 





e2^-Vi_„V(^ 



i-v-)y-2_-^A (13) 



Similarly as r, ii^ is real only for A < 2. Having observed that the energy acquires an imaginary 
part, it is tempting to interpret this as the string breaking. However, this occurs at complex r. 
Our model is therefore stable. 

In contrast to the AdS case HI, the potential in question is written in parametric form given 
by Eqs.® and H13|).^ We can, however, gain some important insights into the problem by 
considering two limiting cases. 

First, let us have a close look at Eq.Q. As noted earlier, the range of A is < A < 2. After a 
short inspection we find that r is a continuously growing function of A. The asymptotic behavior 
near zero is given by ^ 

r = iy^(l-iA(l-V)+0(A^)), (14) 

3 
where p = T^ (^) /(27r)2 . From this it follows that small A's correspond to small values of r. 

The asymptotic behavior near 2 is given by ^ 

-ln(2-A) + 0(l). (15) 

Thus, this region corresponds to large values of r. 

Having understood the correspondence between A and r, we can investigate the properties of 
the interquark interaction at short and long distances. 

We begin with the case of small r. Expanding the right hand side of Ea. (|i;jj) up to the 
quadratic terms in A, we get 

E 



2-Kp 



y|(^l + lA(l-3V)+0(A2)y (16) 



It is unclear to us how to eliminate the parameter A and find _B as a function of r and c. 
^AU the integrals can be evaluated in terms of the beta functions as discussed in the appendix A. 



^The integral is dominated by « ~ 1, where it takes the following form "iJ ^ /„ d«/\/a(l — u) + 6(1 — vY, with 
a = 2(2 — A) and b = A — (2 — A)(3 + 2A). The remaining integral may be found in tables j9|. 



Combining this with (|14() , we find the energy of the configuration as a function of r and c 

E = Q(-^+a,r + 0{r')), (17) 



r 

where Kq = 2^/0~^ and o",) = jcp"^. Thus, we have the expected 1/r behavior at short distances. 
In a similar spirit, we can explore the long distance behavior of E. It follows from H13() that 
in the neighbor of A = 2 the energy behaves as 

E = -^^H2-X) + 0{1). (18) 

Along with the relation (jlSj) , this means that the interquark interaction at long distances is given 
by 



E = g(^jr + 0{l)j (19) 

that is nothing but the desired linear potential. Here we have set a = -^c. 

Having understood the two limiting cases, we can now make a couple of estimates relevant 
to phenomenology. 

It is natural to fix the overall constant q from the slope of the potential at large distances. 
Indeed, the stringy approach is to be most reliable at large distances. From (0) and (fTU]) we have 

Q = —(ca'^y^ ^OM, (20) 

e 

where we have used c « 0.9 GeV^ as it follows from the fits to the slope of the Regge trajectories 



Next, we can estimate the slope of the linear potential at short distances. According to the 

Cornell model (pQ), the slope is the same at all the distances while Eqs. (|19)) and p7j) imply that 
the coefficients in front of the linear terms at large and small distances are different. However, a 
simple estimate of their ratio yields 

— = 87r2er-^(i)«1.24. (21) 

Clearly, the difference in the slopes is not significant for our phenomenological estimates and the 
agreement with the lattice data is very satisfactory at this point. 

Finally, we can compare evaluate the 1/r term in the potential. Phenomenologically, a little 
algebra shows that in (^ the coefficients obey I/ko^ Ri 0.38 GeV^. If we truncate our model by 
keeping only the two terms as in the Cornell model, we find 

"" ^ ■er4(i)c«0.85GeV2. (22) 



Kq IGvr^ 

The value is more than twice bigger than that of the Cornell model. So, this looks disappointing. 
However, the Coulomb-like potential at short distances is controlled by the running coupling 
a sir) and can hardly be predicted within the simplified stringy model we are considering. We 
will come back to discuss this point in the next section. 

To complete the picture, let us present the results of numerical calculations. The parametric 
equation (jUJ predicts a characteristic form for r, as shown in Fig. 2. It has an interesting effect 
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c=0.42 GeV ■ 



c=0.9 GeV 



Figure 2: r as a function of A at c = 0.42, 0.9, and 3 GeV^ 



on the form of the interquark potential in the phenomenologically important interval 0.1 fm < 
r < 1 fm. It is clear that for quite small values of c this interval corresponds to small A's. As a 
result, the approximate formula ()17() is valid. So, if the desired linear behavior holds, it has the 
slope proportional to (Tq. On the other hand, larger values of c result in A's ~ 2. In this case the 
linear behavior has the slope proportional to a. The effect can be seen in Fig. 3 for the window 
0.5 fm < r < 1 fm and the two values of c namely, c = 0.42 and 3 GeV . 
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Figure 3: E/q as a function of r at c = 0.42, 0.9, and 3 GeV^. 



3 Discussions 



In this note we have evaluated the heavy-quark potential using the now standard ideas motivated 
by the gauge/string duality. A key point is the use of the background metric Q which is singled 
out by the observation [7j that it provides, in terms of the same duality, linear Regge trajecto- 
ries. The overall conclusion is that the same background metric results in a phenomenologically 
satisfactory description of the confining potential as well. However, there is a number of open 
problems which are mostly rooted in the heuristic nature of the gauge/string duality in the case 
of pure Yang- Mills theory. In conclusion, we list a few such problems and compare the results 
with theoretical expectations based on more traditional approaches. 

(i) The potential of interest behaves as 1/r at short distances. This behavior is predicted, of 
course, by perturbative QCD as well. On the string theory side, one might fix the overall factor 
g in the world-sheet action ® by fitting the coefficients of the 1/r terms in the energy (|16j) and 



the perturbative calculations. If so, then g ~ ag. By contrast, the AdS/CFT correspondence 
requires g ~ y^ol [H|- Thus, this way of fixing the overall factor does not look satisfactory.^ 

On the other hand, in the case of pure Yang-Mills theory the status of the 1/r term in the 
potential (^ is also unclear. Its numerical value fits the theoretical prediction for the so called 
Liischer term jllj which is derived at large distances. Phenomenologically the fit (^ works at 
small distances as well. Moreover, a straightforward application of the expansion in the running 
coupling as{r) results in a badly divergent series at presently available distances r.^ Thus, it is 
not ruled out that a pure perturbative description sets in at much smaller distances. In view of 
all these theoretical uncertainties, the lack of agreement between the stringy potential (|17)) and 
the Cornell model at short distances might be not so significant. 

(ii) On the lattice, the so called Casimir scaling has been observed ^. The point is that 
if one measures the heavy-quark potential for various color representations of the quarks, at 
large distances the string tension turns out to be proportional to the coefficient of the Coulomb- 
like term. Reproducing the Casimir scaling theoretically is a strong challenge to the model 
building jj^]- Within AdS/QCD, such a scaling is obvious. It is reproduced without fixing any 
parameter. Indeed, the fitting parameter is the overall factor g in the expression for the energy, 
see the discussion above. 

(iii) Power-like corrections to the heavy-quark potential at short distances were studied ear- 
lier mostly in terms of infrared renormalons.^'^ Within this approach, the potential can be 
represented, in somewhat symbolical form, as 

V{r) = (perturbative series) + ^ 6fc r^^A'^^^^ , (23) 

fe=o 

where the coefficients bj. cannot be determined consistently within short-distance physics and 
correspond to the infrared renormalons. In particular, the leading non-perturbative correction 
to the potential is of order [Tl| 

Kon-pert ~ r'^ { CtsC^ )pinst , (24) 

where {agC^) is the gluon condensate and pinst is a typical instanton size. Note that the instanton 
size plays the role of an infrared cutoff. Without such a cutoff, the calculation is divergent. 

In our example, the potential is calculable consistently at short distances. There are no 
uncertainties corresponding to the infrared renormalons. The only exception might be the leading 
renormalon corresponding to /c = in (|23j) . Indeed, by the potential one can understand the 
difference between the total energy and self-energy of the quarks. In the bulk of the paper, we 
were calculating the total energy. Separation of the potential from the self-energy involves then 
the infrared uncertainty in the mass of the heavy quark which corresponds to the renormalon. ^^ 

We close the discussion of the infrared renormalons with a few short comments: 
1. There exists a simple picture of the absence of the renormalons. The key point is that at 
short distances a string doesn't go far away into the fifth dimension or, in other words, it doesn't 
reach a vicinity of the horizon that would correspond to the IR. 



*This difficulty was pointed out in different contexts. See, e.g, |1||10|. 

®For a review, see jj. 
^"For further discussion, see |13|. 
^^At technical level this means that the use of minimal subtraction in Eg. 1121 is not appropriate. 



2. We discuss a class of metric leading to the absence of the renormalons with /c > in the 
appendix B. 

3. Phenomenologically, there is no significant r^ term in the potential. ^^ Thus, the absence of 
the infrared renormalons from the stringy potential can be considered as a success of the model. 

(iv) The AdS/QCD approach provides a natural framework for appearance of a linear piece 
in the potential at short distances. Moreover, as follows from Ea.H21|). the coefficients of the 
linear terms in the short and large distance expansions turn out to be close to each other. 
Phenomenological arguments in favor of such a contribution can be found in J16j . 
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Appendix A 

This appendix collects together some of the formulae that are used in section 2. 
All the integrals can be expressed in terms of integrals of the form 

I{a,b)= [ dvv" {1 - v^Y . (A.l) 

Jo 

Making the variable change w = v'^, we can easily evaluate the integral over w. The final result 

is 

I{a,b) = \B{^,b+l) , (A.2) 

where B denotes the beta function. 

In manipulating the beta functions, we use the relation to the gamma function and the 
formulae [HI 

r(i)=y^, T{l)T{l) = V27r, T{x + l)=xT{x) . (A.3) 

Appendix B 

In this appendix we will investigate the question of whether for a generic form of the warp factor 
h{z) in the metric @ the renormalons with A: > are missing as well.^^ 
Since at z = the metric reduces to that of AdSs, we take 

/i(0) = l. (B.l) 

The Nambu-Goto action and the first integral of equation of motion are given by Eq.© and 
Eq.Q, respectively. From the first integral we get 

^^For further discussion, see |15|. 
^^We assume h{z) > 0. 



where Zq is the inaxhxLum value of z, v = j-, and ho = h{za). 

We compute the energy of the configuration as in section 3. In the process we regularize the 
integral over z by imposing a cutoff e. Finally, we get 



E^ = ^ + E + 0{e), (B.3) 

vre 



where 



^ = ^^ I -1+ / dvv'^ 

TTZo 



Mi-.^f^fr^-i 



(B.4) 



As a result, the potential is written in parametric form given by Eas. HB.2|) and ()B.4|) . 

Assume that we eliminated the parameter Zq and found ii^ as a function of r}^ Now the 
question arises: when does this potential have no terms like r^'^? The answer to this is clear 
from the form of equations. If we transform Zq — > —Zq, then r — > — r for any even function h. 
Meanwhile, the energy transforms as E ^ —E. Thus, E is an odd function of r if /i is even. The 
latter means that in this case the problem of the renormalons with A; > is missing. 
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